Anisotropic Diffusion Limited Aggregation 
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Using stochastic conformal mappings we study the effects of anisotropic perturbations on diffusion 
limited aggregation (DLA) in two dimensions. The harmonic measure of the growth probability for 
DLA can be conformally mapped onto a constant measure on a unit circle. Here we map m preferred 
directions for growth of angular width a to a distribution on the unit circle which is a periodic 
function with m peaks in [— it, n) such that the width a of each peak scales as o ~ l/y/k, where 
k defines the "strength" of anisotropy along any of the m chosen directions. The two parameters 
(m, k) map out a parameter space of perturbations that allows a continuous transition from DLA 
(for m = or k — 0) to m needle-like fingers as k — » oo. We show that at fixed m the effective 
fractal dimension of the clusters D(m, k) obtained from mass-radius scaling decreases with increasing 
k from Ddla — 1.71 to a value bounded from below by D m i n = 3/2. Scaling arguments suggest 
a specific form for the dependence of the fractal dimension D(m, k) on k for large k, form which 
compares favorably with numerical results. 

PACS numbers: 05.45.Df, 61.43.Hv 



I. INTRODUCTION 



Nonequilibrium growth models leading naturally to 
self-organized fractal structures, such as diffusion limited 
aggregation (DLA) yj , have received great interest in the 
recent years not only due to their relevance for various 
physical processes, for example dielectric breakdown 
electrochemical deposition |3j, |J] , and two-fluid Laplacian 
flow but also because such harmonic growth leads 
naturally to one of the most interesting multifractal dis- 
tributions found in nature 



A powerful method for studying such two dimensional 
growth processes is the iterated stochastic conformal 
mapping |g, |9|, LLiJ i which has already been successfully 
applied to generate and analyze DLA |ToL ITlj and Lapla- 
cian |12| growth patterns in two dimensions. This has 
opened the road to address many important questions 
related to pattern formation in DLA, such as the struc- 
ture of the multifractal spectrum of DLA 13], and pro- 
vided the first definite answers for how the hottest tips 
and the coldest fjords grow. Other topics that can be 
investigated using iterated conformal maps include the 
pinning transition in Laplacian growth IjUl . the differ- 
ence between Hele-Shaw flows and DLA 15], as well as 
new topics such as the scaling of fracture surfaces formed 
during quasistatic cracking 16]. 

One of the important questions addressed soon after 
the original discovery of DLA by Witten and Sander [l| 
was that of the effect of the intrinsic anisotropy in lattice 
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models on the shape and fractal dimension of the asymp- 
totic aggregates [13, E, EH H3, W$ ■ For two dimensional 
growth, it was shown that the result of such anisotropy 
in the microscopic attachment probability leads to clus- 
ters which asymptotically have the symmetry of the un- 
derlying lattice (following the argument in Ref E3j this 
actually holds for m < 6, where m represents the coordi- 
nation number of the lattice) , and the fractal dimension 
of the resulting aggregate asymptotically approaches 3/2. 
These results have also been confirmed in recent work 
which used iterated stochastic conformal mapping tech- 
niques to grow the clusters (2^ . 

In the present work we use iterated stochastic con- 
formal mapping techniques to study DLA with m pre- 
ferred directions for growth. Although this naturally 
leads to anisotropic clusters, the present model is fun- 
damentally different from the previous studies on lat- 
tice anisotropy. Our model is rather related to the ex- 
istence of a large scale imposed m-fold symmetry whose 
strength can be tuned. Specifically let us consider the 
case when the harmonic measure for DLA is weighted at 
angle ip between the seed and the location for growth by a 
term W(ip; m, a), where a specifies the angular width of 
the preferred direction. Such a weighting W(tp; m, a) ~ 
exp[— pH.{tp] m, a)] could be due to an imposed exter- 
nal field or to growth on a surface which has an m fold 
symmetry. An example would be dendritic growth in a 
strip [2Gj which can be argued to lie in the m = 1 or 
772 = 2 universality class, the anisotropy increasing as the 
strip is narrowed. 

The organization of the paper is as follows. In Section 
II we describe how we use conformal mapping methods 
together with an angle dependent probability for growth 
P(0; m, k) to study a model corresponding to a real space 
weighting W(ip; m, a). Here is the angle parametrising 
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the unit circle to which the boundary of the growing clus- 
ter is conformally mapped, m is the number of the privi- 
leged directions, and k is an appropriate measure for the 
"strength" of the anisotropy. In Section III we present 
results for the morphology of the resulting patterns as 
a function of m and k, and we derive using scaling ar- 
guments the effective fractal dimension D(m, k) of the 
emerging clusters. We conclude with a discussion of the 
results in Section IV. 



II. MODEL AND THEORETICAL 
BACKGROUND 

In the DLA model proposed by Witten and Sander Q 
the growth of a cluster from a seed placed at the origin 
proceeds by irreversible attachment of random walkers 
released from infinity (in practice, from far away from 
the cluster's boundary). Thus the probability P(s) for 
growth at any point s along the cluster boundary of to- 
tal length L is a harmonic measure and can be written 
P(s) = |(W)(s)|/ J L |(W)(s')Ms', where V(r) is a po- 
tential which outside the cluster obeys Laplace's equa- 
tion \7 2 V = subject to the boundary conditions V = 
on the (evolving) boundary of the cluster and V ~ lnr 
as r — > oo (corresponding to a uniform source of parti- 
cles far away from the cluster). In two dimensions, this 
formulation as a potential problem has been recently ex- 
ploited for studying the time development of DLA based 
on conformal mapping techniques [a ll(l | . 

As discussed in detail in [3, , the basic idea is to 
follow the evolution of the conformal mapping of 
the exterior of the unit circle in a mathematical u-plane 
onto the complement of the cluster of n particles in the 
physical z-plane rather than directly the evolution of the 
cluster's boundary. The equation of motion for $( n )(cj) 
is determined recursively (see Fig.QIa)). With an initial 
condition corresponding to the unit circle in the physical 
plane [oS) = u, the process of adding a new "particle" 
of constant shape and linear scale \Ao to the cluster of 
{n — 1) "particles" at a position s chosen according to 
the harmonic measure is performed using an elementary 
mapping (f>x,e(^>) 
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which conformally maps the unit circle to the unit circle 
with a bump of size VA localized at the angular position 9 
@. The parameter a describes the shape of the elemen- 
tary mapping; following the analysis in |lCj . we have 
used a — 0.66 throughout this paper as we believe the 
large scale asymptotic properties will not be affected by 
the microscopic shape of the added bump. As shown di- 
agrammatically in Fig.^a), the recursive dynamics can 



than be represented as iterations of the elementary bump 
map <fi\ ni g n (Lu), resulting in the convolution representa- 
tion of the conformal map z = & n '(uj) at the nth stage 
of growth as 

$ (n) (w) = <f)\ 1 q 1 o (f) X2fi2 o • ■ ■ o </>A„,e„M . (2) 

where the angle 9 n S [— 7T, tt) at step n is randomly cho- 
sen because the harmonic measure on the real cluster 
translates to a uniform measure on the unit circle in the 
mathematical plane, i.e, 



P(s)ds 



d9 
2?r' 



(3) 



Eq. [3] is crucial to the successful implementation of the 
iterated conformal method as the highly nontrivial har- 
monic measure in the physical plane becomes uniform in 
the mathematical plane. Finally, 
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$(n-l)'( e i0„)|5 



(4) 



is required in order to ensure that the size of the bump 
in the physical z plane is V^o- We note that in the com- 
position Eq. |21 the order of iterations is inverted - the 
last point of the trajectory is the inner argument, there- 
fore the transition from (uj) to <j>( n+1 ) (u>) is achieved 
by composing the n former maps Eq. |21 starting from a 
different point. 
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FIG. 1: (a) Diagrammatic representation of the mappings $ 
and <j>. (b) Change in shape of the probability distribution 
G(£,;m,k) with increasing k for m = 3. 
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Consider the case where the existence of to preferred 
directions in physical space modulates the harmonic mea- 
sure at any point s on the boundary by a probability 
P(ip(s);m,cr) = W{^{s);m,a)/ ^ W{^{s');m,a)ds' . 
Here, ip(s) is the angle parameterisation of the clus- 
ter boundary in the physical space, W{tp] to, a) is 
the modulating weight, and the to- fold periodicity im- 
plies P(il> + 27t/to;to, ct) = P(ip;m,a). The impor- 
tant question is if the weighting W(ip; m, ct) in the real 
space may be represented in the form of a modula- 
tion of the constant measure in the mathematical plane 
Pmath{&) = dO/2-K. Because the angle ip is not invari- 
ant under the conformal map z an answer 
to the question above is not straightforward. Consider- 
ing an ensemble of clusters generated under the influ- 
ence of the same modulation P(-0;to, ct), for each clus- 
ter of n particles ip maps onto a different 6 n (ip), where 
expiif) = <I?(™)(expi6' ra )/|$(™)(expi6> ra )|. It is reasonable 
to assume that averaging over the many patterns above 
an asymptotically (n — > oo) average scale invariant pat- 
tern will appear. For DLA, i.e., in the absence of mod- 
ulation, this pattern is a circle; in the general case, an 
to fold periodic pattern with the same symmetry as the 
modulation is expected to appear. Therefore, we ex- 
pect (#„(?/>)) = f m (ip), where (...) denotes an average 
over clusters, with / m (V0 independent of n and satisfy- 
ing f m (ip + 2ir/m) = f rn (ip)+2Tr/m due to the symmetry 
of the resulting averaged pattern and to the fact that f m 
is an angle. In principle, f m is defined up to an addi- 
tive constant; this can be fixed by specifically requiring 
that at the pattern fingertips j = 0, 1, 2, ...to — 1 one has 
/ m (27rj /to) — 2nj/m. It then follows that the modu- 
lated probability distribution on the unit circle leading 
to such to— fold symmetry patterns obeys 



P, 



till 



(6)=P(f- 1 (e);m,a)df- 1 (0)/de 



(5) 



Thus, we can see that for this case P m ath{9) is itself an 
TO-fold periodic function on the unit circle with peaks at 
the preferred directions 9k — 2?:k/m. 

In this paper, rather than attempting to derive a 
specific P ma th{8) using Eq. we shall directly assume 
such an m-fold periodic measure on the unit circle 
(which based on the arguments above is expected to 
lead to an to— fold symmetry weighting function W in 
the physical space) and study the clusters created us- 
ing the choice P ma th{9) — G(9;m,k)d9, where the pa- 
rameter k ~ is an appropriate measure of the an- 
gular width of the preferred direction in the physical 
plane. The angle-dependent probability distribution on 
the unit circle G(9; m, k) will be normalized such that 
J_ 7r d9G(9;m,k) = 1. Such a distribution biases the 
choice of the location 9, and thus s, where growth oc- 
curs as follows. At step n, the point s for the attempt 
of growth is chosen, as before, based on the harmonic 
measure, i.e., one chooses points 9 n 6 [— 7r,7r) on the 
unit circle with uniform distribution. But growth at s 
is only allowed with a probability G(9 n ; to, k). If the at- 
tempt is rejected, then the previous sequence is repeated 



until a successful trial occurs. We note that obviously 
G(9; m, k) = const corresponds to usual DLA, while an 
explicit dependence on 9 models the existence of privi- 
leged directions. 

Because at this stage we are interested in the general 
features of such a model of anisotropic growth, and not in 
trying to model a specific physical system, we will make 
the simple choice 
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G(0; to, k) = — — cos 
C[k) 

G [ — 7T,7T ) , TO G N , 



where 
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k) 



is the normalization constant (note that it does not de- 
pend on to). 

It is easy to see that G(9; to, k) defined above has all 
the key properties required: for m > it is a periodic 
function of 9 of principal period 2w/m, and thus the num- 
ber to of peaks of G(9; to, k) in [— tt, it) corresponds to the 
number of privileged directions; obviously, both k = in- 
dependent of to and to = independent of k correspond 
to isotropic DLA growth. The exponent k > is a mea- 
sure for the "strength" of selectivity, i.e., the larger k, the 
narrower and higher peaks of G(9; m, k). Therefore, the 
pair (to, k) defines a two dimensional parameter space for 
the analytic study of applied anisotropy to DLA. Though 
the choice given by Eq.[S]is arbitrary, we believe that be- 
cause of universality the key features will be independent 
of the specific form of the function G(9; to, k). 



III. RESULTS AND DISCUSSION 

The model described in Sec. II was simulated as fol- 
lows. The parameter Ao = 10 -3 was fixed because it is 
just setting the microscopic area of an added "particle" . 
For fixed to and k the growth step n proceeds by select- 
ing at random (uniform probability) an angle 9 £ [— 7r, w), 
then comparing G(9] m, k) with a random number r uni- 
formly distributed in [0, 1/C(k)); if r < G(9;m,k), then 
9 n = 9, A„ follows from Eq. and the new map $„ 
follows from Eq. [21 if not, the previous sequence is re- 
peated until a successful trial occurs. All the averages 
mentioned below were done over 100 clusters grown up 
to size N = 20000. An example of clusters with differ- 
ent symmetries (i.e., different values m) grown in this 
manner is shown in Fig.|^a), while Fig.J^b) depicts the 
morphology with increasing anisotropy "strength" (i.e., 
different values k at fixed to). 

It can be seen that the bias introduced by the distri- 
bution G(9; to, k) is indeed producing clusters with the 
corresponding m— fold symmetry and that it is very ef- 
fective: even for small values k, the clusters in Fig. |2^a) 
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FIG. 2: Typical clusters (size N = 20000) grown with (a) 
m = 3, 4, fixed k = 3, and (b) fixed m = 6, but different 
values for k, k = 1, 5, 10, respectively. 



show a clear 3-fold, respectively 4-fold, symmetry. In- 
creasing k (the strength of the anisotropy) leads to a sig- 
nificant reduction in the branched structure of the clus- 
ter, thus in the thickness of the surviving branches, as 
shown in Fig. EJb). Similar results have been obtained 
for all the values 2 < m < 7 and 1 < k < 80 that have 
been tested. 

The results shown in Fig. [21 suggest that the resultant 
patterns have a fractal morphology that depends on m 
and k, and in order to characterize these shapes we will 
focus on the effective fractal dimension D(m, k) obtained 
from the mass-radius scaling. Following the arguments 
in Ref. [HJ, the coefficient = Uf =1 (l + X l ) a in the 
Laurent expansion of $w, 



F[ n) tu- 



F, 



(n) 



■F^lo- 2 - 



(7) 



is a typical length-scale of the cluster; thus, a natural 
choice for the radius of the n-particle cluster is R ~ F 
Assuming that for n >> 1 a scaling law of the form 



(«) 



l/D(m,k) 



(8) 



is found, the effective fractal dimension of the cluster 
can be extracted from a power law fit to the numerical 



data. We note in passing that this scaling law was used 
in Ref. |10| as a very convenient way to measure the 
fractal dimension of the growing DLA cluster. As we 
have anticipated, for all the values m and k the numerical 

(n) 

results for the average coefficient F± show clear power- 
law dependence on the size n, an example being shown 
in Fig. Ota), and the results D(m,k) obtained from the 
power-law fit to the data in the range n > 10 3 are shown 
in Fig. B|». 
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FIG. 3: (a) Average as a function of n for clusters grown 
with k = 1, 10, 40, and 80, respectively, and fixed m = 6 (log- 
log plot). Also shown (dashed lines) are both the limit case 



(n) 



,1/Ddla 



(DLA cluster), where Ddla = 1.71, and 



the proposed lower bound for anisotropic DLA growth [l^| 
F[ n) ~ n V D -min t where D min = 3/2. (b) The effective fractal 
dimension D(m,k), obtained from ~ n 1/D{m - k) , as a 

function of k at fixed m. The points represent the measured 
values, the lines are just a guide to the eye. 

It can be seen that D(m, k) decreases with increasing k 
at fixed m (and with increasing m at fixed fc), and there 
is a certain tendency for saturation at large k. We note 
that, as expected, D(m, 0) ~ -Ddla and that the curves 
D(m,k)aie all above the expected lower limit D m i n = 
3/2 |17|. An exception is the case m — 7 (results not 
shown), where for large k the values D(7, k >> 1) ~ 1.45 
are somewhat below D m i n , but this is most probably 
due to either insufficient statistics (too few clusters), as 
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suggested also by the noisiness of the D(m, k) curves, or 
to the fact that in this particular case the size N = 20000 
is not sufficient to obtain an asymptotic cluster. 

In order to understand these results theoretically we 
will use a simple argument, following Ref. |17|. based 
on the assumptions that (a) for large k the growth of 
the cluster occurs mainly at the tips of the m principal 
branches, and (b) the envelope of the average cluster can 
be approximated by m diamond shaped polygons, like 
the one shown in Fig.^Ja), of opening angles 7 and (3 (in 
general, these angles depend on both m and k) and with 
edges of lengths in the order of R. 
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FIG. 4: (a). Superposition of 10 different clusters of size 
N = 10 5 grown with the same m — 3 and k = 1 but different 
sequences of random numbers. The dotted diamond around 
the arm centered at ij) = shows the approximation for the 
envelope of one arm of the cluster, and the drawing in the 
upper right corner shows schematically the geometry of the 
diamond, (b) Numerical results for ^d-z as a function of the 
scaling variable S = m\f~k. The dashed lines are just a guide 
to the eye for the linear behavior in the range S » 1 and, 
respectively, S > 1. 

Under these assumptions, the rate of growth can be 
written as [l7| 



dN/dR~ iZ*/(3«-«. 



Because the LHS of Eq.EJis LHS ~ RD{m,k)-i^ once the 
angle f3(m, k) is known D(m, k) can be determined from 



D(m, k) = 1 



2tt - P(m, k) 



(10) 



Simple geometry (see the schematic drawing in the top 
right corner of Fig.^Ja)) allows one to write (under the 
assumption that the angles 7(777, k) and /3(m, k) are small 
- which is certainly true for large k and m), 



7 = 



R 



0= 5 -^P=^, (11) 
r r 



On the other hand, the opening angle 7(771, k) is obvi- 
ously fixed by the decay of the probability for growth 
G(9; m, k), and therefore can be estimated as the width of 
the peak of the distribution. Working with the peak cen- 
tered at 9 — 0, assuming large k and small 7, the width at 

half peak probability is given by 1/2 « 1 — (7717/4) 

1 - k (mi/A) 2 . Thus, from Eq.lTTl 



7(777, k) 



i\fk 



P(m, k) 



2(R 



Ci 

lyk ra\fk 



(12) 



where C\ is a constant of 0(1), independent of m, k, or 
the size N of the cluster. Combining Eas.ll0land ll2l one 
thus obtains the following scaling relation for the fractal 
dimension 



D(m, k) ~ 1 7T 
2D(m,k) - 3 = ~C~x 



iVk. 



(13) 



(9) 



Note that it is immediately apparent from Eq. 1131 that 
Dmin = 3/2, while for large k the ratio (D — 1)/(2D — 3) 
should be a linear function of the product S = m\fk 
only, provided S » 1. This prediction can be tested 
against the numerical results. As shown in Fig.^Jb), the 
data collapse is excellent and the function is indeed linear 
when S » 1, confirming our theory. Surprisingly, the 
scaling predicted by Eq. ^| seems to hold down to quite 
small values of k, although these results are beyond the 
scope of our scaling arguments. 



IV. CONCLUSIONS 

Using iterated stochastic conformal maps, we have 
studied the patterns emerging from a model of 
anisotropic (in the sense of the existence of privileged ra- 
dial directions for growth) diffusion limited aggregation 
in two dimensions. In our model, the anisotropy was in- 
troduced via a probability distribution for growth with a 
number 777 of peaks in [— ir, 7r), the width of a peak (giving 
the "strength" of anisotropy) being a tunable parameter 
that allows a continuous transition from isotropic DLA 
growth to anisotropic clusters. We have shown numerical 
evidence that at fixed 777 the effective fractal dimension 
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of the clusters D(m, k) obtained from the mass-radius 
scaling decreases with k from D^la to values bounded 
from below by D m in =3/2. Using simple approximations 
(supported by numerical results) for the envelope of the 
cluster and general scaling arguments, we have derived a 
scaling law involving D(m, k) and successfully tested it 
against numerical results. 

Although the model we have proposed is somewhat 
artificial, it has the advantage that it seems to capture 
most of the general features of an anisotropic growth pro- 
cess while it is still simple enough to allow an analytical 
treatment (to a certain degree). Finally, we note here 
that a system for which the proposed geometry may be 
easily experimentally achieved is the growth of bacterial 
colonies. For such a case, the radial anisotropy can be ex- 
perimentally obtained through the addition of nutrients 



along the privileged directions, and controlled through 
the excess concentration of nutrients along these direc- 
tions in respect to the rest of the substrate. This would 
allow a direct testing of all our numerical and analytical 
predictions. 
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